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I. INTRODUCTION 

Quantum state teleportation is undoubtedly one of the most fascinating developments of quantum information 
processing 

Teleportation of an unknown quantum state is its immaterial transport through a classical channel employing one 
of the most puzzling resources of Quantum Mechanics: entanglement A variety of possible experimental schemes 
have been proposed and few of them partially realized in the discrete variable case involving the polarization state of 
single photons A successful achievement has been then obtained in the continuous variable case of an optical 

field raj . However, the tantalizing problem of extending quantum teleportation at the macroscopic scale still remains 
open. 

Recently, in the perspective of demonstrating and manipulating the quantum properties of bigger and bigger objects 
, it has been shown Q how it is possible to entangle two massive macroscopic oscillators, like movable mirrors, by 
using radiation pressure effects. The creation of such an entanglement at the macroscopic level suggests an avenue 
for achieving teleportation of a continuous variable state of a radiation field onto the vibrational state of a mirror. 

II. THE MODEL 

We consider the situation where an unknown quantum state of a radiation field is prepared by a verifier (Victor) 
and sent to an analyzing station (Alice). Here we shall provide a protocol which enables Alice to teleport the 
unknown quantum state of the radiation onto a collective vibrational degree of freedom of a macroscopic, perfectly 
reflecting, mirror placed at a remote station (Bob) (see Fig. |l|). For simplicity we consider only the motion and 
the elastic deformations of the mirror taking place along the spatial direction x, orthogonal to its reflecting surface. 
Then we consider an intense laser beam impinging on the surface of the mirror, whose radiation pressure realizes 
an optomechanical coupling |j. In fact, the electromagnetic field exerts a force on the mirror proportional to its 
intensity and, at the same time, it is phase-shifted by the mirror displacement from the equilibrium position pC| . In 
the limit of small mirror displacements, and in the interaction picture with respect to the free Hamiltonian of the 
electromagnetic field and the mirror displacement field x(r, t) (r is the coordinate on the mirror surface), one has the 
following Hamiltonian 

H=- J d 2 vP(r,t)x(r,t), (1) 

where P(r, t) is the radiation pressure force P]. All the continuum of electromagnetic modes with positive longitudinal 
wave vector q, transverse wave vector k, and frequency lu — y/ c 2 (k 2 + q 2 ) (c being the light speed in the vacuum) 
contributes to the radiation pressure force. We are adopting the interaction picture with respect to the free Hamil- 
tonian of the electromagnetic field of the continuum and of the field of elastic deformations of the mirror. Following 
Rcf. p|, and considering linearly polarized radiation with the electric field parallel to the mirror surface, we have 



P(r,t)=- 




x {d(k, q)a(k', q') exp[-i(w + uj')t + i(k + k') ■ r] 
+a f (k, g)a f (k', q') exp[i(w + uj')t - i(k + k') • r] 
+a(k, g)a f (k', q') exp[-i(w - ui')t + i(k - k') • r] 
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+a t (k,g)a(k',g')exp[i(u-cy)t- i(k-k') -r]} , (2) 

where d(k, q) are the continuous mode destruction operators having transverse wave vector k and positive longitudinal 
wave vector component g, obeying the commutation relations 

[a(k,q),a(k',q')] = 5(k-k')S(q-q'). (3) 

Furthermore, the electromagnetic wave frequencies ui and a/ are given by u> 2 = c 2 (k 2 + q 2 ) and to' 2 = c 2 (k' 2 + g' 2 ), 
and Ufc, u q denote dimensionless unit vectors parallel to k, q respectively. 

The mirror displacement x(r, t) is generally given by a superposition of many acoustic modes [ jyl ; however, a single 
vibrational mode description can be adopted whenever detection is limited to a frequency bandwidth including a single 
mechanical resonance. In particular, focused light beams are able to excite Gaussian acoustic modes, in which only a 
small portion of the mirror, localized at its center, vibrates. These modes have a small waist w, a large mechanical 
quality factor Q, a small effective mass M |0], and the simplest choice is to choose the fundamental Gaussian mode 
with frequency Q, i.e., 



£(r, t ) = \l^m r m + Pe \ ex p(- r > 2 ) • ( 4 ) 



By inserting Eqs.(0) and (Ef) in Eq.rtll) and integrating over the variable r one obtains 



H = fdkfdqf *> I d q>^L(u k • u fc ,)u, 



8tt 2 V 2Mn _ 
x {a(k, g)a(k', q') exp[-*(w + u')t - (k + k') 2 w 2 /4] 
+a t (k, g)a f (k', g') exp[i(w + uj')t - (k + k') 2 w 2 /4] 
+a(k, g)a f (k', g') exp[-i(w - oj')t - (k - k') 2 w 2 /4] 

+a t (k,g)a(k',g')exp^(w-w')^- (k - k') 2 w 2 /4] } x {be~ int +Pe int \ . (5) 

In common situations, the acoustical waist w is much larger than typical optical wavelengths and therefore we 
can approximate exp { — (k ± k') 2 w 2 /4} w 2 /An ~ <5(k ± k') and then integrate Eq. (|J) over k', obtaining 



H = ~^\hL[dk[dq[dq<^g 



2tt V 2MO _ 
x {— a(k, g)a(— k, g') exp[— i(u/ + 
— a^(k, g)a^(— k, g') exp[i(w + 
+a(k, g)a^(k, g') exp[— i(w — a/)t] 

+a f (k, g)a(k, g') exp[z(w - uj')t}} x {Se~ 40t + SV *} . (6) 

We now make the Rotating Wave Approximation (RWA), that is, we neglect all the terms oscillating in time faster 
than the mechanical frequency f2. This means averaging the Hamiltonian over a time r such that fir 3> 1, yielding 
the following replacements in Eq. (|^) 

exp{±i(u/±w±fi)i} -> — 5{u'±u±n). (7) 

T 

The parameter r is not arbitrary, but its inverse, 1/r = Az^et, is the detection bandwidth, that is, the spectral 
resolution of the detection apparata used at Alice station. 

Since u and to 1 are positive and £1 is much smaller than typical optical frequencies, the two terms S(uj' + u> ± 51) 
give no contribution, while the other two terms can be rewritten as 

—5(w'-w±Sl) = 2irA» det 8tf-q ± )'?^-, (8) 
t c z q± 



where q± — (uj ± fl) 2 /c 2 — k 2 . Integrating over q' we get 
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H = -UAv de 



2MQ 



J dk J dg-J={a(k, q)a f (k, q+) b VUTTl + a(k, g)a f (k, g_ ) 6 f O 



+a t (k,g)a(k,g + )6 t V^Tfi + a t (k, 9 )a(k,g_)6V^ _ f7} , (9) 

where we have used the fact that uj'(q±) =w±fl. 

We now consider the situation where the radiation field incident on the mirror is characterized by an intense, quasi- 
monochromatic, laser field with trasversal wave vector ko, longitudinal wave vector go, cross-sectional area A, and 
power p. Since this component is very intense, it can be treated as classical and one can approximate d(k, q) ~ a(k, q) 
in Eq. ([)]) , where (with an appropriate choice of phases) 



a(k, q )=-id!^^6(k-k )6(q-q ), (10) 

with uj = cyjk 2 + q%. 

Due to the Dirac delta, the only nonvanishing terms in the optomechanical interaction driven by the intense laser 
beam involve only two back-scattered waves, that is, the sidebands of the driving beam at frequencies ojq ± f2, as 
described by 



* = thAvdet \f^u qo \l^i { ^? at (ko ' q+) i + ^? at (ko ' 



g_)6t 



"° + Q a (k 0l q + ) St - (ko , - q _) S I (ii) 



uj q v 



where now q± = \J (u)q ± f2) 2 /c 2 — fcg. The physical process described by this interation Hamiltonian is very similar 
to a stimulated Brillouin scattering [i2fl , even though in this case the Stokes and anti-Stokes component are back- 
scattered by the acoustic waves at reflection, and the optomechanical coupling is provided by the radiation pressure 
and not by the dielectric properties of the mirror. 

In practice, either the driving laser beam and the back-scattered modes are never monochromatic, but have a 
nonzero bandwidth. In general the bandwidth of the back-scattered modes is determined by the bandwidth of the 
driving laser beam and that of the acoustic mode. However, due to its high mechanical quality factor, the spectral 
width of the mechanical resonance is negligible (about 1 Hz) and, in practice, the bandwidth of the two sideband 
modes Av mo d e coincides with that of the incident laser beam. It is then convenient to consider this nonzero bandwidth 
to redefine the bosonic operators of the Stokes and anti-Stokes modes to make them dimensionless, 



"j : ^ ; 'V 27rZ ^"° rfe « ( k o>g-) = 2?r y y"( k o,g-) (12) 



\2-nl\v mode _ Aq _ 

a 2 = 2ttW — a (k , q+) = 2nd —a (k , q+) , (13) 

so that Eq.([ll]) reduces to an effective Hamiltonian 

H eff = -ihxiaxb - a\P) - iU9{a-2p - a|S) , (14) 
where the couplings x an d are given by 



X = qoAu det \ ——J- —J = cob foW „ 15 

2M\l V AVmodeTlOJQ V V 2M\lc z Av mod e 



W^-fl. (18) 



UJQ 



with 00 = arccos(c(7o/^o)i is the angle of incidence of the driving beam. It is possible to verify that with the above 
definitions, the Stokes and anti-Stokes annihilation operators a\ and a 2 satisfy the usual commutation relations 



J i,3 ■ 
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III. SYSTEM DYNAMICS 



Eq. (|14|) contains two interaction terms: the first one, between modes a\ and b, is a parametric-type interaction 
leading to squeezing in phase space |1| , and it is able to generate the EPR-like entangled state which has been used 
in the continuous variable teleportation experiment of Ref. The second interaction term, between modes hi and 
b, is a beam-splitter-type interaction |l5[ ], which may degrade the entanglement between modes a\ and b generated 
by the first term. 

The Hamiltonian ([14]) leads to a system of linear Heisenberg equations, namely 

hi = xS f , (17a) 

b = X a\ - 0a 2 , (17b) 

k 2 =9b. (17c) 

The solutions read 

«i (*) = ^2 [° 2 ~ X 2 cos (0*)] h x (0) + | sin (0*) St (0) - i [ X 6 - X cos (0*)] fit (0) , (18a) 

5(f) = sin (0f) fi|(0) + cos (0f) 6(0) - - sin (0t) fi 2 (0) , (18b) 

- ^2 [X« ~ tfcos (0f)] fit( Q ) _ i sin (et) 6(0) - A [ X 2 - 2 cos (0*)] a 2 (0) , (18c) 



where = yj ' ff 2 — x 2 ■ 

On the other hand, the system dynamics can be easily studied also through the (normally ordered) characteristic 
function ^(/x, v, C), where /z, v, C, are the complex variables corresponding to the operators 6i, b, ai respectively. From 
the Hamiltonian (jl4|) the dynamical equation for 4> results 

$ = X [nv + fj, u — fi — - v*— - u— $ 

\ ov ov* on o/i* J 

C*^-r + C^--^^-^U, (19) 



with the initial condition 

<f>(t = 0) = cxp [-rl|i/| 2 ] , (20) 

corresponding to the vacuum for the modes oi, a 2 and to a thermal state for the mode b. The latter is characterized 
by an average number of excitations n = [coth(?iCl/ 2k bT) — l]/2, T being the equilibrium temperature and ks the 
Boltzmann constant. Then, equation ( |l9| ) has a Gaussian solution of the form 

3> = cxp [-.AH 2 - B\v\ 2 - £\C\ 2 + Cfiv + C\i*v* + TfiC + ^M*C* + 1>K* + T>v*(\ , (21) 

where 

4 2/^2 2 

A(t) = ^ [cos (20*) - 1] - 2^- [cos (0*) - 1] + [1 - cos (20*)] , (22a) 
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B(t) = ^2 - cos ( 2Qi )l + 2 [1 + cos (20i)] ' (22b) 

C(t) = -^3 sin (20*) + ^ sin (0*) + sin (20*) , (22c) 

V(t) = ^ sin (20*) - ^ sin (0*) - sin (20*) , (22d) 

£(t) = |^ [cos (20*) - 1] - 2^f [cos (0*) - 1] + I 1 - cos ( 2et )l - ( 22c ) 

^(*) = f^I t cos ( 26t ) - !] - |s (X 2 + 02 ) t cos " !] + t 1 - cos ( 2Qi )l ■ ( 22f ) 



4 



After an interaction time t, the state of the whole system can be expressed in terms of the normally ordered 
characteristic function as 



p lfc2 = /^ J^J ^*a*, Vj c,t)e-"* l "-""'- |cl "^(-M)A(-^(-o, 



(23) 



where D indicates normally ordered displacement operator. 



IV. TELEPORTATION PROTOCOL 



The idea is to find an experimentally feasible, modified version of the standard protocol for the teleportation 
of continuous quantum variables |l4| , |l5| ], able to minimize the disturbing effects of the beam-splitter- type term in 
Eq. (@). 

First of all, the driving mode is filtered out after reflection on the mirror (see Fig . l it ), allowing only the modes a\ 
and 02 to reach Alice's station. Then Alice performs a heterodyne measurement |l6[ on the mode 02, projecting it 
onto a coherent state \a). Alice and Bob are left with an entangled state for the optical Stokes mode a\ and the 
vibrational mode b, conditioned to this measurement result, i.e., 



Pib = N 



[— ! ^$( M ,,,C,t)e-l-l 2 -IH 2 HCI^ l( _ /i) ^ ( _, )(a |^ 2( _ C) | Q); 

J 7T J 7T J 7T 



(24) 



and the normalization constant is N — (£ + 1) exp [|a| 2 /(£ + 1)] . Denoting with $(/z, v) the normally ordered 
characteristic function associate to the state (p4|), we have 



i>(/x, v, Qe 



^ a * + Ca-\C\ 2 



exp 



A 



T 2 



S + l 



B 



V 2 



£ + 1 



TV 

e + 1 



(fil/ + fJ,*U*) + 



T 



8 + 1 



(a/i — a*/i*) + 



V 



£ + 1 



{av* — a*i>) 



(25) 



Introducing the quadratures 



X ai — 

x b = 



&i + a\ 
b + ti 



Pb = 



ai — a\ 
iy/2 



iV2 



(26a) 
(26b) 



it is possible to evaluate their correlations through Eq.(|25|). In particular, defining v = (X ai , P ai , Xt,, Pb), the 
correlation matrix Tij = (v^Vj + VjVi)/2 results 



r = 



/ A f+i + 2 
e+i 

V 



C + £2. 

^ £+1 



U S+l 



^ £+1 



u £+1 



£+1 




(27) 



B- + i / 

" £+1 T 2 / 



We now employ the standard protocol for the teleportation of continuous quantum variables Jl 4|jl 5| ] . The quantum 
channel between Alice and Bob is established via two- mode entangled state described by the correlation matrix (|27|). 

An input Gaussian state at Alice's side can be fully described by its 2 x 2 covariance matrix T m . Then, the output 
Gaussain state at Bob's side would be characterized by the covariance matrix r ou *. The input-output relation for 
these matrices can be found as follows. In terms of normally ordered characteristic functions wc have 



1 



exp 



- u r°"*u T 



K(u) exp 



1 



in„T 



-uT m u 



(28) 
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where u = (q,p) is the variable vector of the characteristic functions. Instead K is the Fourier transform of the kernel 
in the integral transform mapping the Wigner function of the input state into the Wigner function of the output state 
(see e.g. Ref. 0). In terms of the Wigner function Wab of the state shared by Alice and Bob, it results 

K(u) = J dx A dpAdx B d P B e -^A q -ix Bq +ivAV-iPBP Wab(xa,Pa,x b ,Pb) 



exp 



-4 (<?, -p,q,p)T(q,-p, q,p) J 



Then, it is easy to derive the relations 



(29) 



out 

11 

out 
12 
out 
22 
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= r 

— v l 

— M 



12 



(Tii + 2r 13 
(Ti4 - r 12 -f 
(r 22 — 2r 2 4 



Hr 33 ) 
Hr 44 ) 



r 23 ) , 



(30a) 
(30b) 
(30c) 



Thus, the fidelity of the teleportation protocol can be written, with the help of Eqs.(27) and (BOj) as 



F 



1 + [1 + A(t) + B(t) + 2C(t) -(T- Vfj(£ + 1)] ' 



(31) 



where we have specialized to the case of an input coherent state. In such a case, the upper bound for the fidelity 
achievable with only classical means and no quantum resources is F = 1/2 fisj| . 

The fidelity ( |3l| ) does not depend on the Bob's local operations. In fact these are merely displacements based on 
the Alice's measurement results X + , P_, a, i.e. X b -> X b + \/2X + + V^Reia}^ - T>)/(£ + 1), P b -> P b - v^2P_ + 
V^Imjo!} (J-+T>)/(£ + l). Note that the amount proportional to T j {£ + 1 ) deserves to account for the shifted results 
A + ,P_ obtained by Alice by virtue of the heterodyne detection (see Eq. (p5|)), while the amount proportional to 
T>/(£ + 1) deserves to cancel the diplacement on the b mode caused again by the heterodyne detection (see Eq. (§5b). 

To actuate the phase-space displacement, Bob can use again the radiation pressure force. In fact, if the mirror is 
shined by a bichromatic intense laser field with frequencies vcq and wq + £1, employing again Eq. (to) and the RWA, 
one is left with an effective interaction Hamiltonian 



H act cx be-* + tie 1 *, 



(32) 



where <p is the relative phase between the two frequency components. Any phase space displacement of the mirror 
vibrational mode can be realized by adjusting this relative phase and the intensity of the laser beam. 

Finally, for what concerns the experimental verification of teleportation, that is, the measurement of the final state 
of the acoustic mode, one can consider a second, intense "reading" laser pulse, and exploit again the optomechanical 
interaction given by Eq. (|l4|), where now a\ and a 2 are meter modes. It is in fact possible to perform a heterodyne 



measurement J16| of an appropriate combination of the two back-scattered modes, Z = a% 



if the driving laser 



beam at frequency luo is used as local oscillator and the resulting photocurrcnt is mixed with a signal oscillating at 
the frequency f2. The behaviour of Z(t) as a function of the time duration of the second "measuring" driving beam 
can be derived from Eqs. (|T8|), that is 

Z(t) = oi(t) - &l(t) = i [x + 9} sm(e*)6t( ) 
1 



[e 2 - x 2 cos(et) - x o + x9 cos(et)] Si (o) 

^ [ X 9 + X cos(et) -x 2 -0 2 cos(et)] a\ (0) . 



e 



(33) 



It is easy to see that for cos(9t) = and 0(9 + x) ^S> 9(9 — x) the measured quantity practically coincides with the 
mode oscillation operator 6^(0), thus revealing information on the state of the mechanical oscillator. 



V. RESULTS AND CONCLUSIONS 



Fig. I shows the fidelity @ as a function of the (rescaled) interaction time t for different values of the initial 
mean thermal phonon number of the mirror acoustic mode n. The fidelity F is periodic in the interaction time t 



G 



(see Methods), and we show only one of all possible time windows where F reaches its maximum. The remarkable 
result shown in Fig. ^| is that this maximum value, F max ~ 0.85, is well above the classical bound F — 0.5 and that 
it is surprisingly independent of the initial temperature of the acoustic mode. This is apparently in contrast with 
previous results fl9[| showing that entanglement is no longer useful above one thermal photon (or phonon). This 
effect could be ascribed to quantum interference phenomena, and opens the way for the demonstration of quantum 
teleportation of states of macroscopic systems. However, thermal noise has still important effects so that, in practice, 
any experimental implementation needs an acoustic mode cooled at low temperatures (see however Refs. p0| , p2| for 
effective cooling mechanism of acoustic modes). In fact, we see from Fig. |2j that by increasing n, the useful time 
interval becomes narrower. That means the necessity of designing precise driving laser pulses in order to have a well 
defined interaction time. Furthermore, the time interval within which the classical communication from Alice to Bob, 
and the phase space displacement by Bob have to be made, becomes shorter and shorter with increasing temperature, 
because the vibrational state projected by Alice's Bell measurement heats up in a time of the order of ('J m n) , where 
7 m is the mechanical damping constant. The effects of mechanical damping can be instead neglected during the 
back-scattering process stimulated by the intense laser beam. In fact, mechanical damping rates of about j m ~ I Hz 
are available, and therefore negligible with respect to the typical values of the coupling constants x — @ — 5 X 10 5 Hz, 
and 9 ~ 10 3 Hz, determining the Hamiltonian dynamics (see Methods). Such values are obtained with the following 
choice of parameters: p — 10W, w ~2x 10 15 Hz, SI ~ 5 x 10 8 Hz, Ai/ de t ~ 10 7 Hz, Av mode ~ t^ 1 ~ 10 3 Hz, and 
M ~ 10 -10 Kg, which are those used in Fig. ^. These parameters are slightly different from those of already performed 
optomechanical experiments However, using a thinner silica crystal and considering higher frequency modes, 

the parameters we choose could be obtained. These choices show the difficulties one meets in trying to extend genuine 
quantum effects as teleportation into the macroscopic domain. 

The continuous variable teleportation protocol presented here modifies the standard one of Refs. [ pT|p5| by adding a 
heterodyne measurement on the "spectator" mode a 2 . This additional measurement performed by Alice is important 
because it significantly improves the teleportation protocol. In fact, it is easy to see that if no measurement is 
performed on the anti-Stokes mode, the resulting fidelity for the teleportation of coherent states is always smaller 
with respect to that with the heterodyne measurement. In particular, there is still a maximum value of the fidelity, 
F m ax = 0.80 in this case, independent of temperature, but the useful interaction time interval becomes much narrower 
for increasing temperature. 

It is worth remarking that the present teleportation scheme provides also a very powerful cooling mechanism for 
the acoustic mode. As matter of fact, its effective number of thermal excitations soon after the two homodyne 
measurements at Alice station becomes n e ff = 1+ A + B + 2C — (J 7 — V) 2 /(£ + 1). It reduces to n + 1 in absence 
of entanglement, where I represents the noise introduced by the protocol. Instead, the optomechanical interaction 
for a proper time permits to achieve n e ff = 0.17, i.e., an 80% reduction of thermal noise at once, at the moment 
of Alice's measurement, thanks to the entanglement. To this end, the classical communication and the phase space 
displacement at Bob's site are unnecessary, since they do not affect the state variances. 

In conclusion, we have proposed a simple scheme to teleport an unknown quantum state of a radiation field onto 
a macroscopic, collective vibrational degree of freedom of a massive mirror. The basic resource of entanglement 
is attained by means of the optomechanical coupling provided by the radiation pressure. Here we have shown the 
teleportation of the quantum information contained in an unknown quantum state of a radiation field to a collective 
degree of freedom of a massive object. This scheme could be easily extended in principle to realize a transfer of quantum 
information between two massive objects. In fact Victor could use tomographic reconstruction schemes, again based on 
the ponderomotive interaction (see p3j ) , to "read" the quantum state of a vibrational mode of another mirror and use 
this information to prepare the state of the radiation field to be sent to Alice. The present result could be challenging 
tested with present technology, and opens new perspectives towards the use of quantum mechanics in macroscopic 
world. For example, we recognize possible technological applications such as the preparation of nonclassical states of 
micro-electro-mechanical systems (MEMS) (24|], where the oscillation frequency could be higher and, consequently, 
the working temperature can be raised. 
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Classical channel 



Fig.1 S. Mancini et al. 

FIG. 1. Schematic description of the system. A laser field at frequency too impinges on the mirror oscillating at frequency 
fi. In the reflected field two sideband modes are excited at frequencies u>i — u>o — Q and u>2 — wo + These two modes then 
reach Alice's station. The mode at frequency u>2 is subjected to a heterodyne measurement D2, while the mode at frequency 
uii is mixed in the 50-50 beam splitter BS with the unknown input given by Victor. A Bell-like measurement Dl is then 
performed on this combination and the result, combined with the heterodyne one, is fed-forward to Bob as two bits of classical 
information. Finally, he actuates the displacement in the phase space of the moving mirror. 
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